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Abstract 

The propagation of thermonuclear flame in presupernovae la is 
considered. Front parameters are obtained, some speculations on front 
stability are presented. 



1 Introduction 

The problem of thermonuclear flame propagation in supernovae la still stands. 
Full hydrodynamic simulation requires knowledge of small scale parameters 
of flame: its normal propagation velocity, instability regimes. In this paper 
we show instability manifestations using a toy model. In literature there is 
no clear understanding whether flame front is stable or pulsates under insta- 
bility pp. In the second part we carry out full hydrodynamical simulations 
of flame and obtain flame parameters for the range of densities, a similar 
analysis was performed in [2]. 



2 Toy Model 

Let us consider a simple model for evolution of temperature T and reagent 
fraction c0: 

d t T = Kd 2 x T + u cQ(T-T ), d t c=-u cQ(T-T ), (1) 

where G is a theta-function (a step- function) . The system models deflagra- 
tion burning in solid propellants because two main physical processes that 
drive slow front are presented in it: thermoconductivity and burning itself. 
Medium in supernovae is gaseous, but when flame propagates in the centre 



1 The model was proposed by P.V. Sasorov (ITEP) 
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of the white dwarf, in dense matter (p ~ 10 8 -j- 10 9 g/cm 3 ) density jump is 
low, so hydrodynamical effects are small and evolution matches burning of 
solid medium. Moreover, at Lewis number Le 3> 1, the process of burning in 
supernovae in general is similar to that described to our system. The choice 
of burning rate function is explained below. 

A stationary wave must obey boundary conditions: 

t = 0, x^oo: T = 0, c = l, d x T = 0, t = 0, x -»■ -oo : c = 0. (2) 

The system can be simplified by redefinition of x to put k = 1. We are 
searching for the wave front, so every quantity depends only on £ = x — vt. 
Due to translation invariance we put £ = : T = T (the point of center of 
flame). The system can be easily solved: 

£ > : c = 1, T = T e"<, 

£<0: c = e^/«, T=l- ¥ ^- o e^, (3) 

For more simplification we put v — 1, that means Wo = Tq/{\ — T ). Let us 
finally write down the simplified system and its solution: 

d t T = d 2 x T + u cG(T-T Q ), d t c=-u cG(T-T ), (4) 



£>0: c = 1, T = T e- ? , 

£<0: c = e^, T = l-^e^. (5) 

The stability of such a system under small perturbations can be easily 
considered analytically: 

T = T n . p . + T p , c = c n .p. + c p , (6) 

T p = e*/(0, c p = e^(0- (7) 

After some calculations [3] the following result could be obtained: the system 
is stable when ujq < 6, and perturbations grow exponentially when co > 6. 

Such a system can be easily numerically simulated and full evolution of 
unstable regime could be obtained. The task is set as follows: 

c\t=0 = Ctheor; T\ t=0 = T t h e or, (8) 

where c t h e0 r and T theor are defined in (jSJ) with center at x c and bound condi- 
tions: 

c|x=o = 0, c\ x=L = 1, T\ x=0 = 1, T\ X=L = 0. (9) 
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Table 1: Numerical simulation, v - measured front velocity. 

So we set the exact analytical solutions as initial conditions and watch their 
evolution. The Table [T] shows results of simulations. Solutions could be 
split into two groups: "flame" and "therm", "flame" - is the evolution 
as stationary flame front with constant velocity, this regime exists when 
ujq < 6 (according to the table of results), "therm" describes flame decay, 
like evolution under thermoconductivity without burning (example of such 
evolution is shown in Fig. [I]). Zero burning rate when T < T forbids flame 
ignition after "therm" regime, so it clearly cuts only unstable evolution. 
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Figure 1: Front positions for Uq = 8 at different time moments: to < t% < 
t 2 < h. 
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The analytical predictions are in a very good agreement with numerical 
simulation of the model. So this model and its modifications may be used 
for theoretical study of unstable flame fronts. 

3 One-dimensional flame properties 

Let us consider full hydrodynamical evolutions of flame in presupernova la. 
Our goal is to study flame acceleration and deflagration to detonation transi- 
tion [1] . Typical parameters of medium in the center of WD are the following 
[5]: T ~ 10 9 K, p ~ 10 9 g/cm 3 , chemical composition - 12 C+ 16 0. Let us sup- 
pose that only 12 C remains in chemical composition. For given conditions the 
following relations hold: Pr <C 1, Le ^> 1. It means that thermoconductivity 
is the only diffusion mechanism that matters in this case. The coefficient of 
thermoconductivity is the sum of two parts: electron conductivity [6] and 
radiative conductivity [7j. We consider the only one nuclear reaction for 
approximation (also this approach gives ability to study physical effects): 

12 C + 12 C ^ 24 Mg*. (10) 

Here Mg* means the excitation state, it decays through 3 channels: with 
n, with p, with a; and we take into account the sum of rates, which could 
be found in [8]. Its caloricity is q\ = 5.5 ■ 10 17 erg/g. The reaction is the 
first in the network, it occurs between two highly charged nuclei (it is one of 
the slowest), so we could suppose that it determines the whole rection rate. 
Futher burning could be indroduced in our model by changing caloricity. For 
burning up to 56 Ni it will be q2 = 9.2 • 10 17 erg/g. Nuclear screening should 
be taken into account because T = {E cou \)/kT ~ 1. 

The problem is formulated in the following way: a full one-dimensional 
hydrodynamical system of equations with thermoconductivity and nuclear 
reactions is considered. It is solved by our numerical code FRONT (where 
an implicit numerical scheme with a Newton iterations solver is used j9]). 
Initial region of calculation is filled uniformly with 12 C at given T and po- 
Right wall should be free for stream. Left wall is heated by the linear law 
T = T + t(T\ — Tq)/t. Where r obeys r 3> L/c s (L - the size of region 
of interest, c s - sound speed). Such conditions lead to deflagration wave 
ignition by the hot wall. The sequental flame positions are shown on Fig. 
[2j Table [2] shows results of normal flame speed determination for different 
initial density p . 

It should be emphasized here that all velocities obtained are correct only 
for one reaction in network 12 C+ 12 C. Use of full nuclear network changes the 
speed radically [TO] . 
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Table 2: Measured deflagration flame front parameters for reaction 12 C+ 12 C. 
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Figure 2: Sequental temperature distributions (for different physical time). 



4 Conslusions 

The toy model was presented for several puproses: first it clearly shows how 
instability of flame front manifests: it leads to front destruction, second, the 
model could be used in theoretical speculations and in academic studies. Full 
hydrodynamic simulations shows stable front propagations at all considered 
densities (with no evolution shown by "therm" regime in toy model), so we 
could state that thermonuclear flame front is stable in all density range from 
2 • 10 s g/cm 3 to 2 • 10 9 g/cm 3 , the reason for that should be explored. Flame 
front paratemeters were obtained in hydrodynamical simulations. Front ve- 
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locities differs greatly from determined in [2]. The reason is in very simple 
nuclear network: only one reaction is taken into account (for detailes see |10j). 
The work is supported partly by Federal Program "Scientific and pedagogical 
specialists of innovation Russia" contract number 02.740.11.0250, partly by 
the Russian Foundation for Basic Research grant RFBR 10-02-00249-a, by 
SNSF SCOPES project No. IZ73Z0-128180/1, and "Dynasty" foundation. 
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